Investigating the impact of a large river and its surrounding contextual conditions on pedestrians’ summer thermal perceptions in a Cfa-climate city

Thermal comfort studies are usually employed to find subjective thermal responses [indicated by neutral temperature (NT), i.e. the temperature with no thermal stress] of residents from a region towards thermal environments. According to the recently published works in the literature, NTs are affected by many factors, such as geographical location and microenvironments. To elucidate the origins of these effects, the impact of microenvironment elements around a water surface on pedestrians’ thermal perceptions was systematically investigated in this work. The Fujiang River (FJR) in Mianyang City was taken as the sample site. The municipal meteorology station is located next to the site by around 2.5 km. By performing meteorology measurements combining questionnaires, it was found that the riverside NT (indicated by physiologically equivalent temperature, PET) of Mianyang in the summer of 2023 was 21.4 °C. The relationship between the distance from the water (DFW) and NT was quadratic linear. The same phenomenon took place by using either PET or Universal Thermal Climate Index (UTCI) indexes. Meanwhile, the meteorological contexts also affected NTs, including relative humidity (RH) and air velocity (Va). Regarding RH, the NPET increased from 15.2 °C (RH = 50%) to 26.9 °C (RH = 90%). In contrast, the NPET dropped from 23.0 to − 50.6 °C when the Va increased from 0.2 to 2.5 m/s, respectively. From our analysis, it was demonstrated that human thermal responses are significantly affected by both the microenvironmental and meteorological backgrounds around the water surface. Our work provides valuable insights for the proper use of water surfaces in urban design for adjusting thermal comfort.

The global warming and urban heat island (UHI) problems undoubtedly affect people's life quality and city livability 1 .In recent years, extremely high temperatures have frequently emerged in many areas around the world.This effect has caused a series of social-related issues including poor health 2 , thermal discomforts 3 , and psychological stress 4 .Consequently, residents need to rely on mechanical equipment for thermal adjustment 5 , which causes a large amount of energy consumption 6 , resulting in energy crisis 7 and economic pressure 8 .Therefore, improving thermal comforts in energy-free ways has been an integral task of the sustainable development process 9 .Various strategies relating to this process have been deployed in recent years.
Several approaches have been explored to reduce the dissipated heat stress.The studies related to this effect in outdoor spaces are defined as outdoor thermal comfort (OTC 10 ).More specifically, the goal of OTC studies involving subjects is to improve the comfort level of humans.Subjective responses towards certain thermal conditions are exported, which is supported by professional techniques.This approach involves subjective factors and objective (thermal) environments.Subjective perceptions are usually numerically modelled.This can be evaluated by the thermal sensation vote (TSV) factor, which typically classifies thermal levels into 7 scales ranging from • A large water surface (river) was found that might be microclimatic influential in the city as the sample.
• A variety of points around the waterbody were selected that would be meteorologically affected by the dif- ferent degrees of the field measurement.
• The data at each site by field surveys were collected (meteorological parameters by sensors, subjective infor- mation by questionnaire, and physical environments by WinSCANOPY 28 ).• The data of different aspects were associated with each other to find their statistical associations.

Methodology The sample city
The Cfa refers to the climate type with the subtropical humid meteorology 29 .It is cold in winter and hot in summer with high RH in warm seasons.This climate is very common in southern China 30 .Mianyang is located in the southwest of China is a typical Cfa city.With the continuous urban environment degradation, such as pollution 31 and natural land shrinkage 32 , extreme weather conditions have frequently appeared in recent years 33 .The summer of Mianyang city spanned from June to Early September every year; The annual hottest month has a maximum temperature above 35.0°C (Fig. 1  www.nature.com/scientificreports/year 2022 has experienced the hottest summer of Mianyang in history (Fig. 1).Therefore, avoiding the impact of extreme weather is currently emerging as an urgent task.
Mianyang is located in the north of Sichuan, China (Fig. 2), and is geographically surrounded by mountains.A big river (the Fujiang River, FJR) goes through the city.The width of the river section in the study area is 228.0 m, the maximum depth is above 26.0m, and it flows from northwest to southeast (Fig. 3) 35 .The cooling effect of the water surfaces has been widely confirmed to be effective 36 .The cooling intensities are strongly dependent on the local physical properties 37 and could be also affected by the attributes of FJR.FJR supplements the local water demands while adjusting the city climate.In the present work, the sites around the riverbank were selected as samples for field measurement.According to the literature, it has been confirmed that waterbodies and site geometry affect the thermal comfort levels.As a result, the sites located closer to the water surfaces 38 and/or with lower openness 39 are cooler.
Therefore, points with different distances from the water surface (DFW) on the bank were selected for performing the measurements.There were a total of 16 points chosen, which were divided into two groups (Groups A & B).The test site was around 2.5 km from Mianyang Meteorological Station (Fig. 3).Group A refers to sites canopied by trees, whereas the B were mostly open.All points in each group ranged on a line being vertical to the river.They were defined as Points 1 (within 5 m of the river) to 8. The neighboring points were 10.0 m to each other.
The representative Groups A (canopied) and B (open) along the FJR were selected as samples.According to the degree of vegetation cover, the sites can be classified as open (without canopy cover) or closed (densely  www.nature.com/scientificreports/covered by trees) 41 .Their coverage was indicated by the sky view factor (SVF 42 ).The SVF is a parameter developed by WinSCANOPY 28 supporting the forestry and agriculture studies.It defines the site openness by the visibility of the sky in a hemisphere, ranging from 0 (fully covered) to 1 (totally open).The values of this factor could be acquired from the software in WinSCANOPY.The points in Group A have lower SVF than B; the lands on two lines were covered lawns.They were defined as A1 to A8 and B1 to B8. Figure 4 presents the images of the physical characteristics of all points (locations, surrounding environments, & SVF).www.nature.com/scientificreports/

Meteorology measurement
The measurements were carried out between 16 and 23 July 2023, which was the annual hottest time.All measurement days were cloudless.There were two sets of Testo 480 supporting the meteorology data collection including T a , V a , RH, and globe temperature (T g ).They were implemented by different sensors and their properties are listed in Table 1.Every point was measured by one full daytime (8:00 a.m.-8:00 p.m.).The two points were fully measured each day.The measuring information of each point is presented in Table 2. Table 3 lists the analyses of standard deviation for the meteorology data.

Questionnaire
Subjective information of all interviewees was obtained by questionnaire, which was carried out simultaneously with meteorology measurement.The questionnaire was divided into two parts according to the research content (Table 4).The first part was about subjective perceptions.TSV, spanning as − 3, cold; to + 3 hot; and + 4 very hot) 44 .The − 4 (very cold) was ignored as nobody would vote that in such hot weather.The setting followed ISO 44 and ASHRAE 45 .Volunteers were suggested to vote with a resolution of 0.1 in TSV 27 .The second part was about the personal information that would affect thermal perceptions including clothing, gender, height, etc.It took about 5 min for the subjects to understand and to answer all questions.For a large number of questionnaire sheets collected, some volunteers who frequently visited the riverbank were repeatedly interviewed.A total of

Thermal comfort indices
The thermal comfort is affected by various parameters.To address this issue, complex indices have been proposed in the literature.They contained both meteorological parameters and volunteers' individual factors (i.e., clothing isolation).Normally, PET 14 and UTCI 46 have been frequently used.PET is based on the Munich Personal Energy Balance Model (MEMI), which simulates the thermal conditions of the human body in a physiologically relevant way.It can be defined as the T a at which the human body's thermal budget in a typical indoor environment (without wind and solar radiation) is in balance with the skin temperature under complex outdoor conditions 14 (MRT = T a , work metabolism = 80 W, clo = 0.9, etc.).UTCI can be defined in another s (V a = 0.5 m/s, work metabolism = 135 W, RH = 50% etc.) 46 .
A large number of works in the literature have demonstrated that the perceived comfort values of PET and UTCI are very similar 16 .The common equation for MRT can be expressed as follows 47 .
The signs in the formula refer to MRT, T g , T a , V a (v); whereas D and ε are the diameter of the globe (D = 0 0.15 m) and its absorption rate (0.95), respectively.The values of the indices have been calculated by RayMan Model 14 (Fig. 5) and developed by Matzarakis et al. 48.Either of the three indices has been used in OTC studies.UTCI was developed in 2009 by international co-operation between leading experts in the areas of human thermophysiology, physiological modelling, meteorology, and climatology 49 .By comparing the neutral temperature and neutral range of the PET and UTCI indices, significant differences were found between the neutral temperature and the neutral range in summer and winter 50 .The increase of MRT by 10.0 °C contributed to a TSV ascent of 0.9 51 .

Data analyses
Multiple linear regression (MLR) is usually utilized for data analyses in OTC studies 52 .It is an analytical method that creates a linear relationship between a dependent variable and one or several independent variables.The MLR was used in this work to explore the linear correlation between thermal perception and thermal environment under the action of different microenvironmental factors.The impact of riverside microenvironmental factors on people's thermal perceptions was explored here, which is affected by the NT variation.NTs of all points needed to be calculated by MLR.The TSV and PET were available as dependent and independent variables.

Ethical approval
This study was ethically approved by Southwest University of Science and Technology (23zx7107).

Data description
Figure 6 describes the daily weather of all measurement days (8:00-20:00).As can be seen, they were totally very hot (above 25.0 °C on average).Generally, it got cool from 16 (30.0°C) to 20 (27.0 °C) in July, yet getting warm afterwards, peaking at 31.0 °C (23 July.).Figure 7 shows the personal information of the volunteers when they received the questionnaire.
Figure 8 illustrates the daily PET fluctuation of all 16 points.Each image contains 4 points measured within 2 days.Generally, Points A had lower values than the performances of tree canopies.A1 was cooler than B1 by nearly 10.0 °C at 11:00.As the daily weather changed, all PET values were irregularly changed.Most of them peaked around 17:00 (e.g., 50.0 °C for B2 & 57.0 °C for B3).Nevertheless, the points measured on the same day expressed similar trends.For instance, B6 increased from 25.0 °C (8:00) to 46.0 °C (12:00), and then reduced till 14:00 (40.0 °C).A6 expressed a similar trend despite the different extreme values.Three characteristics for all MRT = (T g + 273.15) 4  PET can be derived: (1) the existence of higher values at open points; (2) two points measured on the same day exhibited similar waves; (3) various changing regulars were recorded at all measurement days.Waving regulars of UTCI are expressed in Fig. 9. Their trends were generally similar to PET.For example, A1 waved similarly with B1, since measured on the same day (increasing from 8:00 to 12:00, bottoming at 14:00, peaking again at nearly 16:00).Open sites were warmer than canopied ones as well, such as B3 being hotter than A3 by 15.0 °C at 17:00.The values of UTCI at the points of all days fluctuated differently.
Figure 10 depicts the Boxplots of TSV ranges at all points.Points measured on the same day were placed nearby for comparison.In fact, the two nearby points had close scopes, which resulted from the effects of the daily weather conditions.As can be seen from the image, Points A had lower values in total, as far as canopies  (2.9) on average.Similarly, they ranged differently, while remarkable scopes emerged at Points 5 and 6, nearly from -1 to 4. There were few regulars about their variation revealed.The subject responses were affected by current meteorology and physical environments commonly.Thus, it is hard to explain the properties of the TSV scopes by the site features.

General neutral points
The neutral temperature The subjective thermal responses were associated with the two OTC indices in separation.As can be seen from the images (Fig. 11), TSV increased as getting warm.The increase of PET by 10 °C caused a TSV rise of around 1.3; in contrast, that of UTCI by the same value increased TSV by 1.9.They were both significantly correlating (R 2 = 0.6 for PET & 0.6 for UTCI).The NTs obtained from the images were 21.4 °C (PET) and 26.1 °C (UTCI).
As shown in Fig. 12, the TSV of Lines A (shaded) and B (open) are linearly correlated with PET and UTCI.The NT values indicated by PET (NPET) of shaded space and open space were 22.3 °C and 22.7 °C, respectively.In contrast, the NUTCI values (NT indicated by UTCI) were 27.7 °C and 26.2 °C, respectively.For the NPET, the correlation of Line B (R 2 = 0.7) was greater than that of Line A (R 2 = 0.5).The sheltered space had a slightly lower neutral temperature (NT) than the open space.For NUTCI, the correlation of Line B (R 2 = 0.6) was greater than that of Line A (R 2 = 0.5).The sheltered spaces had higher NTs than the open spaces.

Neutral points for wind and relative humidity
Figure 13 displays the relationship between the RH and humidity sensation vote (HSV) and the relationship between V a and wind sensation vote (WSV).A linear regression of RH against HSV was performed.From our analysis, it was found that the HSV gradually went from very dry (HSV = − 3) to very wet (HSV = 3) as RH increased.At y-axis = 0, the RH was 45.2%, indicating that at this point in time one feels that the RH was neither dry nor wet (HSV = − 3).V a was linearly regressed against WSV.It was found that as the V a increased, the WSV also gradually went from no wind (WSV = 0) to high wind (WSV = 4).However, the correlations were poor in both two models (R 2 < 0.1 & 0.1).

The neutral temperature in various contexts
The neutral temperature varying for distance from the river In this work, it was assumed that the water surfaces, as a microenvironmental factor affecting meteorology, would affect subjective thermal perceptions (represented by NTs).The intensities of the influencing varied for distance.Table 6 shows NPET and NUTCI at the 16 points.NTs increased with DFW.The highest value of NPET was found in A4, A5, and A8 (NPET = 26.5 °C), and the lowest value was found in A1 (NPET = 18.8 °C).The highest NUTCI value appeared in A8 (NUTCI = 31.1 °C), and the lowest NUTCI value appeared in B1 (NUTCI = 24.3°C).Therefore, it can be inferred that people had lower NTs if being closer to the river.Nevertheless, the river would www.nature.com/scientificreports/have limited this perceptive effect since its affecting intensity was steady when DFW was above 40 m (Fig. 14).Similar trends were witnessed for NUTCI (Fig. 15).

The neutral temperature varying for relative humidity
According to the above-mentioned analysis, it was found that people responded differently towards DFW.As far as it concerns the cooling properties of water surfaces (air humidity adjustment 53 ), the finding might result from the variation of RH.This study assumed that NTs would vary for RH.Therefore, the associations between TSV and PET were analyzed by various RH contexts.The RH values were divided into several groups regarding values in measurements (35-95%, 10% inter value).A variety of NT values towards certain RH ranges were acquired, while some unusual results were obtained (NTs below 0).This might result from special or extreme meteorological conditions.The unusual data were excluded in the further analyses, outputting Fig. 16.The overall trends showed insignificant variation.NTs increased from 15.0 to 27.0 °C as RH raised from 50.0 to 90.0%.

The neutral temperature varying for air velocity
People's thermal perceptions would be affected by other meteorology parameters.NTs were also clustered into various scopes regarding the wind conditions.The measured V a ranged from 0.02 to 2.95 m/s.They were grouped with an interval of 0.3 m/s.With some extreme values excluded (NT values below 0, same as Fig. 17), Fig. 17 was derived.People were increasingly thermal tolerated as getting windier, which could be explained by the lower NT values.The NTs decreased from 23.0 (PET) and 27.0 (UTCI) °C at V a of 0.1 m/s to 15.0 (PET) and 14.0 (UTCI) °C at V a of 1.7 m/s.Significant linear correlations were also found.

Discussion
The impact of various physical (trees and a water surface) and meteorology factors (RH and V a ) on subjective thermal perception in a Cfa climate city (as assessed by NTs) was systematically explored here.The NTs (indicated by PET & UTCI) varied for all points and showed a polynomial correlation with contextual parameters.Similar phenomena were witnessed in meteorological context variation.This resulted from various reasons.People's thermal responses varied for site openness.As far as PET is concerned, the sites covered by trees had lower NTs.The impact of SVF on OTE/C has been widely examined in the literature.Sites more densely canopied were usually cooler.Song 54 found that MRT increased by 2.7 °C linearly since the SVF rise by every 0.1 point.Pedestrians might adapt to locations with various SVF thermally.However, NTs in this study were non-linearly correlating with SVF.This implies that the impact of openness on subjective perceptions was complex, which is in line with the reported results of Xiong et al. 27 .and Zhang et al. 55 and Zhang et al. 23 .They have found, subjective thermal perceptions varied for SVF, yet not always linearly correlating.
People's thermal perceptions varied also for DFW.The water surface undeniably affects the meteorology conditions for surrounding areas 38 , depending on their distances 25 .More specifically, the sites located too far from the water surfaces would be beyond the affecting ranges of them.Du et al. 56 found that land surface temperature kept generally steadily as moving farther from the river although it linearly rose at river nearby areas.Consequently, the heat index would polynomially rise with the DFW increase.The polynomial changes of subjective thermal responses (NTs) towards DFW were revealed here, which highlight the effect of water surface cooling and individual thermal adaptation (Fig. 13).
The thermal responses of humans varied for meteorology.This study has found that NTs varied for RH and V a .The NTs reduced for V a , which was confirmed by Zhang and Lin 57 and Hou 58 .Generally, windy conditions would cause less heat stress.The wind would be improve the water evaporation 59 , which is cooling effective 60 , reducing heat stress.In other words, people would feel cool under warm if being windy 61 .This phenomenon meets the results of the present study.Furthermore, their NTs varied differently in RH levels.The RH has been proven to be effective on examining thermal comforts 62 .Spaces with higher RH would have denser water molecular weight, which absorbs heat energy in the atmosphere, and lead to a reduced heat stress 63 .In other words, personal heat stress would be reduced by the water molecules during hot seasons 64 .Nevertheless, the underlying origins of this effect are still in controversy.There should be positively or negatively influential 65 .
Overall thermal responses (NTs) were found in this work.They were indicated by PET and UTCI.OTC in Mianyang and some cities with different climates.have been frequently studied in recent years and compared in Table 7.The lowest and highest NTs were found by Cheng et al. 66 and Xiong et al. 27 , respectively.A relatively low NT value was found in this work.This might result from the water cooling effects and human thermal adaption.According to the literature, the test points mainly focused on the representative area of crowd activity 67 or the point with the largest flow of people 66 .On this basis, the study expanded the selection range of distance from the water (DFW ≤ 80 m).The changes in the neutral temperature (NT) and distance from water (DFW) were studied, and the range of influence of the water on the heat perception of humans on the shore was determined.
Some shortcomings exist despite the study's important findings.Only polynomial correlations between NTs and contextual factors were found.In fact, the involved factor (DFW) was perceptively thermal effective as its original cooling effects 38 .Other influential factors, e.g., SVF and parameters indicating individual trees, were insignificantly explored.Correspondingly, the influential distances would be affected by river properties, such as width and depth, which were not considered.These factors might be perceptively influential as they were broadly proved meteorological impactive (SVF; tree crown dimension; river width).Nevertheless, some unusual phenomena were reported in Sections 3.3.2and 3.3.3.The NTs were found to be below 0 at some special

Figure 1 .
Figure 1.Histogram of extreme temperature values recorded in Mianyang in the latest decade in summer.

Figure 2 .
Figure 2. Locations of the city and the studying areas 35 (processed by Adobe Photoshop 40 ).

Figure 4 .
Figure 4. Locations, surrounding environments, and sky view factor images of all points (processed by Adobe Photoshop 40 ).

Figure 5 .
Figure 5. RayMan Model for the calculations of PET and UTCI.

Figure 6 .
Figure 6.Weather changes of all measurement days by mean T a from the weather station.

Figure 7 .
Figure 7. Personal information about volunteers.

Figure 8 .
Figure 8. PET fluctuations for all points during the measurement periods.

Figure 9 .
Figure 9. UTCI fluctuations for all points during the measurement periods.

Figure 12 .
Figure 12.Linear correlation between TSV and thermal indices on the two lines.

Figure 16 .
Figure 16.Linear correlation between NTs (indicated by UTCI & PET) at different RH conditions.

Table 1 .
43operties of all sensors of Testo 48043.

Table 2 .
The measurement date and time.

Table 3 .
The standard deviation of meteorology for the two groups.

Table 4 .
Questionnaire survey.original data generated or analysed during this study are included in this published article as the supplementary (the excel file named as 'The original data 1214') file.

Table 5 .
Number of questionnaires issued and collected.